******************
**** Figure 4 ****
******************

* Run table b2 to get estimates W *
do "$dofiles/table b2" 

* Plot
coefplot W, keep(treatmentalldummy2 treatmentalldummy3 treatmentalldummy4 treatmentalldummy5 treatmentalldummy6 treatmentalldummy7) vertical yline(0, lcolor(black)) recast(bar) barwidth(0.9) lcolor(gs0) fcolor(*.5) color(gs6) ///
ciopts(recast(rcap)) citop  coeflabel(treatmentalldummy2="P-Bas" treatmentalldummy3="N-Bas" treatmentalldummy4="P-Req" treatmentalldummy5="P-Attr" treatmentalldummy6="N-Req" treatmentalldummy7="N-Attr") format(%9.2f) ///
msymbol(t) levels(95) graphregion(color(white)) addplot(scatter @b @at, ms(i) mlabel(@b) mlabpos(2) mlabcolor(black)) ///
ylabel(0  0.1  0.20  0.30  0.40  , angle(0)) ///
ytitle(ITT Effect on Collection Efficiency, color(black)) title(Treatment Effects During Intervention, color(black)) subtitle("June 2017 - August 2017")  ///
groups(treatmentalldummy2 & treatmentalldummy3  = "{bf:BASIC}" treatmentalldummy4 & treatmentalldummy5 = "{bf:PIDA}" treatmentalldummy6 & treatmentalldummy7 = "{bf:NIDA}")


graph export "$filepath/output_logs/figure 4.pdf", replace
